Correlation-induced metal insulator transition in a two-channel fermion-boson model 



00 

o 
o 

s : 

<: 
■ 

<n : 



o3 



i 

c 

o 
o 



On 

O 

(N 
O 
oo 
O 



G. Wellein, 1 H. Fehske, 2 A. Alvermann, 2 and D. M. Edwards 3 

J Regionales Rechenzentrum Erlangen, Umversitat Erlangen- Nurnberg, 91058 Erlangen, Germany 

2 Institut fur Physik, Ernst- Moritz-Arndt-Universitdt Greifswald, 17489 Greifswald, Germany 

3 Department of Mathematics, Imperial College London, London SW7 2AZ, United Kingdom 

We investigate charge transport within some background medium by means of an effective lattice 
model with a novel form of fermion-boson coupling. The bosons describe fluctuations of a correlated 
background. By analyzing groundstate and spectral properties of this transport model, we show 
how a metal-insulator quantum phase transition can occur for the half-filled band case. We discuss 
the evolution of a mass-asymmetric band structure in the insulating phase and establish connections 
to the Mott and Peierls transition scenarios. 



The way a material evolves from a metallic to an in- 
sulating state is one of the most fundamental problems 
in solid state physics. Besides band structure and dis- 
order effects, electron-electron and electron-phonon in- 
teractions are the driving forces behind metal-insulator 
transitions (MITs) in the majority of cases. While the 
so-called Mott-Hubbard MIT |l[ is caused by strong 
Coulomb correlations, the Peierls MIT [2] is triggered 
by the coupling to vibrational excitations of the crystal. 
Both scenarios are known to compete in subtle ways. 

The MIT problem can be addressed by the investi- 
gation of generic Hamiltonians for interacting electrons 
and phonons like the Holstein 01, 0|, Hubbard [3, [j| 
and (quarter-filled) t-J models |6(, or combinations of 
these [7|. These models have been proved to describe 
MIT phenomena for the half-filled band case, in partic- 
ular for one-dimensional (ID) systems known to be sus- 
ceptible to the formation of insulating spin-density-wave 
(SDW) or charge-density-wave (CDW) broken-symmetry 
groundstates [8|. On the metallic side of the MIT, 
charge transport then takes place within a "background 
medium" that exhibits strong correlations, which antic- 
ipate the SDW or CDW on the insulating side. In that 
particle, as it moves, creates local distortions of 
substantial energy in the background. These distortions 
may be parametrized as bosons. They are able to relax, 
with a rate that depends on the system properties but 
also on the proximity to the MIT. 

In order to model such a situation the authors recently 
proposed a simplified transport Hamiltonian J^, [lfj| 



H = H h 



(1) 



where Hi, = —tbYlu cjcj(bj + bj) describes the boson- 
affected nearest neighbor (NN) hopping of fermionic par- 
ticles (cj) 11]. In (JTJ) the particle creates a boson {b\) 
on the site it leaves and destroys a boson on the site it 
enters. Thereby it generates a "string" of local bosonic 
fluctuations with energy ujq ■ Cutting the string, the 
A term allows a boson to decay spontaneously. A unitary 
transformation 6j t— > bi + X/loq eliminates the boson re- 
laxation term in favor of a free-particle hopping channel, 



H 



f = ~tf E<i,j) 4 C i Wlt11 l f 

the original one. As a result 



with tt = 2Aib/wo, in addition to 



H ^ H = H b + H f + coo b i bi 



(2) 



and the physics of our model is governed by two pa- 
rameter ratios: the relative strengths of the two trans- 
port channels (tf/tb) and the rate of bosonic fluctua- 
tions (wo/tb) -1 . The model has been solved numerically 
in the one-particle sector and revealed — despite its seem- 
ing simplicity — a surprisingly rich "phase diagram" with 
regimes of quasi-free, correlation and fluctuation domi- 
nated transport [l(| • In this case the spinless Hamilto- 
nian |2]) covers basic features of the more complicated t- 
J, Hubbard or Holstein models in the low doping/density 
regimes, but is much easier to evaluate. 

Whether our two-channel transport model likewise de- 
scribes a quantum phase transition from a metallic to an 
insulating phase at certain commensurate band fillings 
remained an important but open question. The free hop- 
ping channel Hf will clearly act against any correlation 
induced charge ordering that might open a gap at the 
Fermi energy Ep. Strong bosonic fluctuations, i.e. small 
uj , will also tend to destroy CDW order. On the other 
hand, a tendency towards CDW formation at half-filling 
is expected for large ivo/tb by perturbative arguments, 
yielding an effective Hamiltonian with nearest-neighbor 
fcrmion repulsion. In some respects this is evocative 
of the quantum phase transition in the spinless fermion 
Holstein model, which for large phonon frequencies and 
strong couplings can be mapped on the XXZ model 
The XXZ model undergoes a Kosterlitz-Thouless transi- 
tion at the spin isotropy point. 

Further evidence of a quantum phase transition comes 
from the investigation of simplified versions of |TJ: (i) a 
coarse-grained model with only one bosonic oscillator for 
the ID infinite system, and (ii) a ID model based on a 2- 
site cluster. The exact solution of (i) and an approximate 
solution of (ii) both exhibit quantum phase transitions 
at A = between two Fermi liquids, with different Fermi 
surfaces, non-Fcrmi-liquid behavior persisting down to 
zero temperature at the critical point in case (ii) 
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The MIT is a subtle quantum mechanical problem, 
however, that requires non-approximative investigation 
schemes. Therefore, in this work, we apply unbiased 
numerical techniques to study the competition between 
itinerancy, correlations and fluctuations for the ID half- 
filled band case, using the full model |T|) without any re- 
strictions. To this end we employ exact diagonalization 
in combination with kernel polynomial expansion meth- 
ods, adapted for coupled fermion-boson systems [HI]. 
The computational requirements are determined by the 
Hilbert space dimension Dh = (jv )( iv b )' wnere N ^ s 
the number of lattice sites, N e counts the fermions, and 
Nf, is the maximum number of bosons retained. Typically 
we deal with Dh of about 10 11 . 

Let us start with the discussion of the photoemis- 
sion (PE) spectra. The spectral density of single- 
particle excitations associated with the injection of an 
electron with wave- vector k, A + (k,u>) (inverse PE), 
and the corresponding quantity for the emission of 
an electron, A~(k,iv) (PE), are given by A ± (k,uj) = 



Here c^T 



and 



\ipo) is the groundstate in the iVe-particle sector while 
\4>n) denote the n-th excited states in the N e ± 1-particle 
sectors with excitation energies = — Eq. 

Figure [T] displays the wave-vector resolved single- 
particle spectra in the regime where distortions of the 
background are energy-intensive, i.e. the boson fre- 
quency loq is high. If the free transport channel is dom- 
inant (tf = 5 - upper graph), the occupied (unoccu- 
pied) band states, probed by PE (inverse PE), give rise 
to an almost particle-hole symmetric absorption spec- 



trum A+ (uj - E F ) ~ A K ^(E F 



Thereby the main 



spectral weight resides in the uppermost (lowest) peaks 
of A~^ (A~k) in each if-sector. The corresponding "co- 
herent" band structure roughly follows the — 2tfCosK 
tight-binding band. Satellites with less spectral weight 
occur near the Brillouin zone boundary predominantly, 
as a result of mixed electron-boson excitations with total 
wave vector K. At T — the Fermi energy is obtained 
from J2k J~L A K(u)duj = N e = N/2 (half-filling, no 
spin). We see that there is no gap between A± v ^ 2 and 
A± x / 2 a t Bp. Moreover the spectral weight of both peaks 
is almost one, i.e. a particle injected (removed) with 
K = Kp = ±7r/2 propagates unaffected by bosonic fluc- 
tuations. The system behaves as an unusual metal. 

If we decrease A (t//if, ratio) at fixed u>o, we enter 
the regime where boson-assisted transport becomes im- 
portant (see lower graph of Fig. Q]). At about A c (wo = 
2) ~ 0.1 a gap opens at K = ±tt/2 in the PE spec- 
tra. The gap increases as A < A c gets smaller, but its 
magnitude shows no finite-size dependence (to demon- 
strate this we included the N = 16, iV& = 9 data for 
K = ±7r/2). Most notably Ep lies inside the gap re- 
gion, signalising the transition to the insulating state. 
The MIT is correlation induced. Since A is small, distor- 
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FIG. 1: Photoemission (black) and inverse photoemission 
spectra (red) for the half-filled band case with loq — 2 at 
tf = 5 (upper panels) and tf = 0.01 (lower panels), where 
N — 12, Nb = 15. Dashed lines give the integrated spec- 
tral weights, e.g. S^(uj — Ef) = f£ duj' A%^{oj' — Ef), where 
Sk = Sx(— oo) + S~k(oq) — 1, and ~}2 K Sk = N. Here and in 
what follows periodic boundary conditions were used, leading 
to discrete K(N) wave numbers. All energies are measured 
in units of tb = 1, and u is rescaled with respect to Ef- 



tions of the background cannot relax easily. Accordingly 
the band structure is strongly renormalized. We observe 
that now A 



k \ 



Ep) ~ A^_ k (lu — Ep) and expect a 
perfect doubling of the Brillouin zone for N — > oo. Map- 
ping our model ((T|) for u ^> tf , t h to a (XXZ-like) spin 
model, S z — ► — S z symmetry is broken, reflecting the ob- 
served broken particle-hole symmetry: The highest oc- 
cupied states belong to an extremely flat quasiparticle 
band, whereas the lowest unoccupied states are much 
more dispersive [bH ]. 

The correlated band structure can be understood by 
"doping" a perfect CDW state (Fig. [2]). To restore the 
CDW order a doped hole can be transferred by a coher- 
ent 6-step process of order 0(t|/wg), | o • •) — ► | * o • ) — » 
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FIG. 2: Doping a perfect CDW, states with one particle re- 
moved (left panel) are connected by a 6-step hopping process 
(see text), whereas a 2-step process (right panel) relates states 
with an additional particle. 



| * * o ) — > | * =f< ) — ^ | o * *) — > | • ■ o), where in steps 
1-3, three bosons (*) are excited, which are consumed in 



steps 4-6 afterwards [lfj, [15( . In this process the fermion 



(o) becomes correlated with the background fluctuations. 
We note that such a coherent hopping process, in which 
the particle propagates and restores the background, ex- 
ists even for the case A = where transport is fully 
boson-assisted. In contrast an additional electron can 
move by a two-step process of order O(t%/ujo). Conse- 
quently the electron band is much less renormalized than 
the hole band, and the mass enhancement is by a fac- 
tor 0((£;,/wo) 4 ) smaller. Note that the mass-asymmetric 
band structure that evolves here is correlation induced. 

That the observed MIT is indeed correlation induced is 
corroborated by the weakening and Anally closing of the 
excitation gap if the boson energy luq is reduced at flxed A 
(see Fig. [3]). In this way the ability of the background to 
relax is enhanced, fluctuations overcome correlations and 
the system turns back to a metallic state. At the same 
time the spectral weight is transferred from the coherent 
to the incoherent part of the (inverse) PE spectra, espe- 
cially for K away from Kp = ir/2 where the lineshape is 
affected by rather broad bosonic signatures. 

The CDW structure of the insulating state be- 
comes apparent by investigating the particle-particle 
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FIG. 3: (Inverse) Photoemission for u>o = 0.5 and A = 0.01, 
i.e. tf(\,Lo ) = 0.04. Again N = 12, but now N b = 15. 
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FIG. 4: Particle-particle (% e e) respectively particle-boson 
(Xeb) correlation functions [upper graph], and weight of the 
m-boson state (|c m | 2 ) [lower panel] in the groundstate of the 
half-filled ID two-channel fermion-boson model ([5}. 



hJ2i( n i n i+j) and particle-boson Xeb(j) 



Xee(j) 

W~ ^2i( n i"i+j^i+j) correlation functions, where rii = ctq. 

In Figure|4]the even-odd modulation of the charge den- 
sity away from a singled out site i of the first particle is 
clearly visible. We note that the charge structure factor, 
Sc(tt) = ^ 1) J (0|K- 1/2) { nj - 1/2) |0), increases 

by a factor of about two in going from A = 0.1 to A = 0.01 
for lo q = 2 (cf. Tab.© 0. 

In the CDW, where e.g. the even sites are occupied, 
every hop of a fermion excites a boson at an even site. 
This gives a large contribution to Xeb(j) at even sites in 
addition to NN sites \j\ = 1 (middle panel). Since the 
CDW involves only few bosons (see lower panel), this 
is the dominant contribution in first order of tb/u>o-> an d 
explains why the boson density is large at sites with large 
fermion density, although the hopping term tj, creates 
bosons at the neighboring sites of a fermion. 

The charge oscillations become rapidly suppressed by 
increasing A, but there is still a reduced charge density 
at the particle's neighboring sites, which enhances the 
mobility of the carrier. Accordingly the boson density is 
enlarged (suppressed) at the NN sites (site) of the parti- 
cle. Clearly Xeb(j) is small V|j| if ujo ~> h, tf because of 
the high energy cost. As expected the fluctuation domi- 
nated regime is characterized by a large number of boson 
quanta in the groundstate. The position of the maximum 
in the boson weight function |c m | 2 is shifted to slightly 
larger values as A increases, i.e. the correlations weaken. 

Further information on the groundstate properties can 
be obtained from the kinetic energy parts Ekinj/b — 
(0|-ff//b|0). Table U shows that boson-assisted hopping 
becomes the major transport mechanism at small A. 
On the metallic side of the MIT the creation and an- 



4 



A 




*/ 




Ehin,b 


D 


S c (ty) 


5.0 


0.5 


20.0 


-149.190 


-13.915 


80.336 


0.0417 


0.1 


0.5 


0.4 


-2.745 


-14.878 


4.054 


0.0419 


5.0 


2.0 


5.0 


-37.040 


-5.291 


20.290 


0.0417 


1.0 


2.0 


1.0 


-7.008 


-6.214 


5.146 


0.0425 


0.1 


2.0 


0.1 


-0.568 


-7.036 




0.0561 


0.01 


2.0 


0.01 


-0.023 


-7.125 




0.1056 



TABLE I: Groundstate properties of the two-channel trans- 
port Hamiltonian © at half-filling (N = 12, N b = 15). 



that the CDW groundstate is a few-boson state. Obvi- 
ously the system is more susceptible to CDW-formation 
at large boson frequency u> (small transfer amplitude tf) , 
keeping the boson relaxation A = tftOo/2th fixed. This is 
the limit of an effective fcrmionic system with (instanta- 
neous) Coulomb repulsion. Recall that as a consequence 
of the correlation-induced CDW state, a band structure 
with a very narrow valence and broad conduction band 
evolves, different in nature from simple two-band models. 
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nihilation of bosons opens a coherent transport channel 
in the regime where strong correlations persist in the 
background. The Drude weight D, obtained from the 
f-sum rule —D = \{E kin j + E kin , b ) + J ° c ' a reg (uj)dui, 
serves as a measure for this coherent transport. Here 
OYe 9 M = J2 n> o ~ u n) is the regular part of 

the optical conductivity (with current j = jf + %). 

At the MIT point D vanishes for the infinite system. 
At the same time the optical gap opens. In the insulating 
phase the optical response is dominated by multi-boson 
emission and absorption processes. Thus the spectral 
weight contained in the regular part of <r(u>) is enhanced. 

A small boson frequency allows for large fluctuations 
in the background, i.e. many bosons in our model. This 
supports transport via the 6-step process on the one hand 
but, as in the one-particle sector [10|, also limits the mo- 
bility of a particle by many scattering events. Neverthe- 
less D is expected to stay finite even for uj /tb — > 0. 

To summarize, the two-channel transport Hamiltonian, 
introduced for studying the dynamics of charge carri- 
ers in a correlated/fluctuating medium, has previously 
only been properly analyzed for a single carrier [l(J. In 
this limit the model may capture some of the physics of 
2D high-T c superconducting cuprates [17[ or 3D colossal 
magnetoresistive manganites [18| . Here we focused on 
the metal insulator transition problem at finite particle 
density, in particular in ID at half-filling, which might 
be of importance, e.g., for the ID CDW MX chain com- 
pounds [8| . Since in this case the problem is of the same 
complexity as for the quarter-filled t-J z -J± or spinless 
fermion Holstein models we make use of elaborate nu- 
merical techniques in order to avoid uncontrolled approx- 
imations. From our finite-cluster study we have strong 
evidence that the model exhibits a quantum phase tran- 
sition from a metallic to an insulating state. The MIT is 
driven by correlations, like the Mott-Hubbard transition, 
but in our case true long-range order develops because a 
CDW state is formed. This might point towards a Peierls 
transition scenario. The Peierls instability, however, is 
most pronounced in the adiabatic limit of small phonon 
frequencies, with many phonons involved in establishing 
the CDW (lattice dimerization) . By contrast we find 
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